Mode coupling and evolution in broken-symmetry plasmas. 
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Abstract 

The control of nonlinear processes and possible transitions to chaos in systems of interacting 
particles is a fundamental physical problem. We propose a new nonuniform solid-state plasma 
system, produced by the optical injection of current in two-dimensional semiconductor structures, 
where this control can be achieved. Due to an injected current, the system symmetry is initially 
broken. The subsequent nonequilibrium dynamics is governed by the spatially varying long-range 
Coulomb forces and electron-hole collisions. As a result, inhomogeneities in the charge and velocity 
distributions should develop rapidly, and lead to previously unexpected experimental consequences. 
We suggest that the system eventually evolves into a behavior similar to chaos. 
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Plasmas are of interest in subjects as diverse as astrophysics and the design of quantum 
solid-state nanostructure devices .— t2»M They exhibit a variety of nonlinear phenomena, even 
close to equilibrium, including instabilities and chaotic processes on different scales.-^ The 
development of strong turbulence, characterized by Porkolab and Chang 5 as a "stochastic 
collection of nonlinear eigenmodes" , is a general, and still puzzling, feature of plasmas. The 
Coulomb interaction between carriers plays the crucial role in producing such a collection of 
coupled modes. Due to the very complex dynamics, the ability to control the coupling and 
evolution of nonlinear eigenmodes is a challenging problem. 

Recent progress in optical phase control allows the production of plasmas in semiconduc- 
tors with a well-controlled charge density and, more importantly, a well-controlled current 
densityM^^ The control of the initial current density is achieved by the quantum interference 
of a one-photon transition (light frequency 2u, with the field phase fau) and a two-photon 
transition (light frequency u, with the field phase across the fundamental band gap 
E g . At nonzero A0 = (p2 U — 20^ the symmetry of the injected distribution in momentum 
space is broken, and a macroscopic current with a speed Uq = t> e |sinA0| is injected in a 
direction parallel to the sample surface. The maximum speed of the injected electrons, v e , 
is determined by to and E g , reaching 10 3 km/s for 2hu — E g about 100 meV. 

Studies of nonequilibrium electron processes in semiconductors^ show that the entire 
dynamics is complex even for a uniform electron density. When current is injected, the re- 
sulting separation of electrons and holes leads to strongly nonuniform Coulomb forces. Here 
we consider situations where these forces determine, rather than just perturb, the devel- 
opment of the charge and current density patterns that can lead to possible nonlinearities 
and instabilities. The system we study theoretically is a multiple quantum well (MQW) 
structure, consisting of up to tens GaAs/Al x Gai_a;As periods, each of thickness on the order 
of 15 to 30 nm, grown along the z direction. At photon energies where carriers are injected 
only in the GaAs layers, the total thickness w of the region, that is the number of periods 
multplied by the period width, where the plasma is produced in typical MQWs can be on 
the order of 0.1 fim, still considerably less than the spot size of the exciting laser beams 
and the light absorprion length; the fact that allows to treat all single quantum wells as 
equivalent electrostatically coupled layers, neglecting the direct motion of elecrons between 
the wells. The injected carrier densities are typically Gaussian in the two-dimensional coor- 
dinate r = (x,y), given by N e ^ h (r,t = 0) = N exp (— r 2 /2A 2 ), (e for electrons, h for holes) 
where A is the spot size, and No is the maximum total injected two-dimensional density for 
all quantum wells, which is proportional to the total number of single quantum wells and 
can be on the order of 10 13 cm -2 . Nq is the concentration parameter in our analysis. As 
a result, the three-dimensional density distribution can be modelled^ as uniform along the 
z-axis, with: 



For this reason we treat the density and velocity distributions in the (xy) plane only. 

The in-plane electric field depends on an integral over the charge density —eN c (r,t), 
where N c (r, t) = N e (r, t) — Nh (r, t) , e is the fundamental charge, and is given by 
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where the model Coulomb kernel J^c (d) = d/ (d 2 + w^) 3 ^ 2 takes into account the width 
of the system and simplifies the calculations by avoiding the singularity at d — 0. Here 
e is the background dielectric constant. The parameter wc is on the order of structure 
width, where for wc -C A the results are not sensitive to the kernel behavior at d C A. 
The field E(r, t) is very sensitive to the details of the carrier dynamics, since even relatively 
small changes in N e (r,t) can strongly modify it. For example, even if N e ^{r,t) are taken 
to be slightly separated identical Gaussian profiles, Eq.([2]) shows that E(r,t) is strongly 
nonuniform. Nonuniformities in the field and in the velocities and the charge patterns 
mutually enhance each other. This process is our main interest. 

To study the nonlinear dynamics, we employ a hydrodynamic model for the dynamics of 
injected charge currents and densities, and include the possibility that an external electric 
field is also present. In hydrodynamic models one avoids requiring the details of distribution 
functions by constructing approximate, closed sets of equations involving conserved and 
slowly-varying quantities such as charge, momentum, and energy densities. In the effective 
mass approximation, closed equations in the range of parameters we consider can be obtained 
for the velocity and density— For simplicity, we assume that the holes in the injected plasma 
are not moving, which does not qualitatively influence our results^ due to a small effective 
mass ratio of electrons and holes. The injection typically occurs on a time scale of 50- 
100 fs. We take this as instantaneous, and treat it as preparing our initial state. Since the 
timescales of interest are much shorter than electron-hole recombination times, the dynamics 
is governed by the continuity equation for the electron density and the Euler equation: 

^ + V(iV e u) = 0, 

^ + ( uV ) u + — = e ( E + jj) _ _ £ ( 3 ) 

dt m e N m e r ch N r e ' 

where E is a time-dependent external electric field. Here and below the r and t— arguments 
are omitted for brevity; P is the pressure, m e is the electron effective mass, the weakly 
concentration-dependent r e h describes momentum-conserving drag 13 due to the Coulomb 
forces during electron- hole collisions, 14 r e is the relaxation time due to external factors, such 
as phonons^ and disorder. Here we consider the effect of this drag only, assuming r e ^> r h 
for a clean sample and electron energies too low for intense phonon emission.— The electron- 
hole drag and the Coulomb forces, being coordinate-dependent, increase the inhomogeneity 
in the charge density. 

To obtain the solution of equations we use a finite basis set, following the Galerkin 

method, and convert Eqs.Q to a system of ordinary differential equations. The expansion 
has the form: 

^rnax ^max 

N = E Nm*n, ut = E + U t (t) , (4) 

n n 

where i — x, y is the Cartesian index. To improve the convergence, we include known 
functions of time Ut (t) in the right- hand- side of Eq.(jH) for velocities. These functions can 
be obtained by solving the equations of motion in the rigid-spot approximation 16 - where 
the electron puddle moves with uniform velocity u = (U x (t), U y {t)) while keeping its initial 
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Gaussian shape. The initial distribution A^ exp (— r 2 /2A 2 ) suggests the eigenstates of a 
harmonic oscillator ^n(x, y) = i[> ni (x)il> n2 (y) as the basis set of the expansions with: 

i> n (x) = 1 H n (x/A) e~ x2 / 2A \ (5) 



where H n (x/A) is the Hermite polynomial of the nth order, and the double-index n = 
(rii,n 2 ). The basis functions satisfy the conditions for norm and derivatives: 



ipn 2 (x)ip ni (x)dx = A8 nitTl2 , (6) 
V2Ai/j' n (x) = y/mj} n -i(x) - Vn + ltp n+1 (x). 

In this basis, the matrix elements for the components of the Coulomb integrals C^.^ (T2J) are 
given by: 

= j~ 2 J J *m W K§ (r - r) * K (?) d *rd\. (7) 
Taking into account (EI), the equations of motion can be written in the operator form: 



dNw 
dt 



-L_lj2 N n (ulP 2 L 12 P x + v^P x L X2 P 2 ) - U (i\ - h) N m 



dv x p 2 u- N h II N h ( pE dU\ 

% = E^.*-£-^p 2 --^ + ^ — , (8) 

dt em e V m ' k ^ r ch N r ch N \m e dt J 

where the equation for du^ /dt is similar to the latter. Here we assume E parallel to the 
x-axis, and put U y = and U = U x for the current injected along the x-axis. The small 
terms (uV) u and VP/(m e N) in the Euler equation have been neglected; the justification 
of this approximation will be given later in the text. We have put: 

/OO (^r£ 
j\) n% (x)i/; ki (x) (x) -j-, 

f°° . . . dx f°° , , , dx 
Im = / iPmi (x) -r / ^m 2 0*0 -r- ( 9 ) 



A 7-oo r '" 2 v ' A 

The operator L 12 = q + £2 ~~ ^1 ~~ ^2, where the ladder operators t v and ^ act on the 
corresponding index, for example: i 2 Pi = \fkiP ni ^ l -\, mi - For the problem we consider here, 
the initial conditions are: A^(0) = \^N 5 ni fi5 n2 fi and w^(0) = 0, where N£ is nonzero only 
if n\ = n 2 = and remains constant in time. Some of the interesting gross quantities that 
can be calculated with these equations will be analyzed below. 

The electron-hole drag makes % dependent on all components of the velocity. Despite 
this complication, Eqs.([8]) can be solved directly in the case of vanishing long-range Coulomb 
forces. The resulting charge density has the form: 

N c (T,t) _U tX ex {_ ± e -rV2AA (1Q) 



iV e (r,0) A 2 *V r eh 
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FIG. 1: (Color online.) Inhomogeneity parameters of the charge density (main plot) and velocity 
(inset) pattern for two regimes of the Coulomb forces. The parameter p = Q, p \T e h is shown near 
the plots, p = 16 corresponds to the extremely weak damping, while for p = 0.25 the damping is 
relatively strong. Here T p \ = 0.9 ps. The functions presented in the plots are universal in the sense 
they do not depend on the initial speed of the puddle Uo. 

The inclusion of long-range Coulomb forces leads to a much more complex dynamics. In 
Eqs.® the Coulomb matrices C— couple a given velocity component to all density 
components N% allowed by symmetry. In turn, the density evolution depends on all possible 
products of components of velocity and density. Therefore, a perturbation in one component 
can cause a growing response in a large number of them. The temporal behavior of the 
system is determined by three independent time scales: drag-induced r e h; the plasma period 
T p i = 27r/O p i, where Q p \ is the two-dimensional plasma frequency for the Gaussian density 
distribution,— 1 ^ = ir 3 ^ 2 N e 2 /(4em e A); and the timescale of the external E(£). We use 
the parameter p = f2 p iT e h to characterize the relative effects of the long-range Coulomb forces 
and drag. 

In our simulations we use GaAs m e = 0.067m, where m is the free electron mass and the 
dielectric constant e = 12; we take iVo = 10 13 cm -2 and A = 1 /iin. These parameters result 
in a plasma period T p \ close to 0.9 ps, which is considerably larger than the injection time. 
The parameter wc in the Coulomb kernel is taken as 0.1 A. The basis set includes 32 states 
for each coordinate, giving a convergence^ in the time interval of interest < t < T p \/2. 
We consider different values of p in the experimentally achievable range. 13 

To trace the evolution in the inhomogeneity of the charge density and velocity patterns, 
we study £ c and £ u , defined to be ratios of gross quantities: 




(11) 



5 



{u x -Ufd 2 r = J (Vu x fd 2 r, 

that serve as characteristic lengths. Taking into account that the spatial inhomogeneity 
(internode distance) of the function ip n (x) scales at large n as n _1//2 , the number of harmonics 
forming the corresponding pattern scales as A 2 /£ 2 (t) or A 2 /£ 2 (t) if the distributions are 
strongly nonuniform. As one can see in Fig.l, both patterns, especially the density, become 
strongly inhomogeneous and the role and the number of the higher harmonics grows with 
time. Therefore, we expect that the spatial scales of the variations in the density and velocity 
will rapidly decrease. Eventually, a hydrodynamic description will fail, as stochastic behavior 
develops."^ 




FIG. 2: (Color online.) Patterns (in arbitrary, same for both columns, units) of charge density 
N c (x,y) (upper row) and velocity u x (x,y) (lower row) at t = T p \/2. Left column: p = 0.25, right 
column: p = 16. The density has N c (x, —y) = —N c (x,y) symmetry. In the upper row, larger bow 
at x > corresponds to N c (x,y) > 0. The velocity satisfies the condition u x (—x,y) = u x (x,y). 
For u x (x,y) maximum values are achieved at the wings \x\/A close to 3, y = 0. Minimum values 
are achieved at \x\/A close to 1, y = 0. 

The underlying charge density is shown in the upper panel of Fig. 2 where we plot the 
profiles N c (x,y,T p i/2). The lower panel shows the velocity u x (x, y, T p \/2). The profiles 
have a rather complex form, showing that the distributions of both quantities are strongly 
inhomogeneous. We calculate the mean spot displacement: 

^)=Ett h*n(x,y)d\ (12) 

n t 

where N t = ttNqA 2 is the total number of injected electrons. The displacement x(t) has a 
complex time-dependence, after initially evolving simply as U t. Even at later times x{t) is 
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FIG. 3: (Color online.) (a) Mean spot displacements for free spot propagation. Dashed line 
corresponds to the linear undamped oscillations. The displacement of the spot is on the order 
of C/oTpi/4, that is 20 nm for typical Uq = 100 km/s and assumed here T p \ = 0.9 ps. (b) Mean 
spot displacements for the forced oscillations. Dashed lines correspond to the linear oscillator in 
Eq. (|13p . The frequency coefficients are k = 4 (off-resonance) and k = 1.2 (close to the resonance). 

proportional to Uq if all other parameters are kept the same. We show in Fig. 3 the mean 
displacement x(t) denned in Eq. ffl2|) for two different cases presented in Fig.l: considerably 
(p = 0.25) and very weakly (p = 16) damped regimes. An astonishing result is the absence 
of the plasma oscillations even close to the clean limit with p = 16. On the timescale 
of half of the expected oscillation period T pl , the spot becomes strongly inhomogeneous 
with harmonics up to ni,U2 < 20 contributing to the results. Therefore, no well-defined 
oscillations occur. In all cases considered, the maximum of x(t) ~ Uq min(T p i, r e h) is much 
less than A, and therefore the Vu and VP-originated terms in the Euler equation can be 
neglected. 

As another example of this unusual behavior, we present the results for the clean system 
(p = 16) driven by an external field E(t) = E sin(kQ p \t) for the same initial Gaussian density 
distribution as above, but with no current injection (Uq = 0). Here the inhomogeneity 
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develops more slowly than if current were injected, since x(t) increases as t 3 rather than as 
t at the initial stage of the process. Nonetheless, the x(t) is considerably different from the 
expected for a linear oscillator: 



with xq = eEo/rriefl^, due to the fact that the excitation of the higher Hermite- Gaussian 
modes strongly influences the response to the external field, as shown in Fig. ([3b). For 
a system driven close to resonance (k = 1.2), the difference between the full and linear 
oscillator behavior is less than for k — 4, since near resonance the uniform external force is 
more important than the interactions. 

To conclude, the macroscopic dynamics of optically injected currents in clean semiconduc- 
tor multiple quantum wells is strongly inhomogeneous and nonlinear, due to the nonuniform 
long-range Coulomb forces that develop. These forces arise following the initial breaking of 
the symmetry by the injected electron puddle velocity Uo, which leads to a separation of 
electrons and holes that produces the nonuniform macroscopic Coulomb interaction. Due 
to the coupling of the Hermite-Gaussian modes through conservation of charge, the charge 
density becomes nonuniform on progressively smaller spatial scales. In contrast to what 
might be expected, it does not show well-defined plasma oscillations. The complex charge 
and current density patterns develop on a time scale on the order of a quarter of the plasma 
oscillation period characteristic of the given carrier density and puddle size. The length 
scales characterizing the spatial inhomogeneities in density and velocity decrease rapidly, 
and, in the terminology of Porkolab and Chang a turbulence regime will likely develop. 
These systems will provide a new laboratory example of plasmas with controlled non-linear 
behavior, and likely a transition to a stochastic regime. 
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